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ABSTRACT 

We present the results of Ho- monitoring of the BL Lac object OJ 287 with the VUT during seven epochs in 2005-08. We were able 
to detect five previously undetected narrow emission lines, /1/16548, 6583[NI1], A6563lla and AA6716, 6731 [Sll] during at least one 
of the epochs and a broad Ho- feature during two epochs. The broad Ho- luminosity was a factor ~ 10 lower in 2005-08 than in 1984 
when the line was previously detected and a factor ~ 10 lower than what is observed in quasars and Seyfert galaxies at the same 
redshift. The data are consistent with no change in the position or luminostity of the Ha line in 2005-08. The width of the Ha line 
was 4200 ± 500 km s"' , consistent with the width in 1984. 
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The BL Lac object OJ 287 (z - 0.306) has received some 
attention after the discovery of the recurring outbursts ev- 
ery ~ 12 years in the optical data over the past century 
dSillanpaa et al.lll988l) . This observation has led many authors to 
consider a binary black hole model for OJ 287 (Sillanpaa et al. 
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Valtonenll200 7r The exact mechanism of the outburst differs 
from one model to another. Some models attribute the flux in- 
crease to c hanges of the direction of the relativistic jet (e.g. 
lKatzlll997l) . other models to the impact of the secondary black 
hole onto the accretion disk of the primary followed by t idal ef- 
fects and increased luminosity of the j et (e.g.lL ehto & V altonenI 
11996) . The most detailed model so far dValtonen 2007) suggests 
a very high mass for the primary black hole M - 1.8 x 10'" Mq 
and a very eccentric orbit for the secondary (e - 0.66). To ex- 
plain the aperiodicity of the outbursts the precession of the orbit 
has to be very strong, A(p = 40° per orbit, in this model. 

Testing these models was so far concentrated mainly on the 
timi ng of the outbur sts with broadband photometric observations 
(e.g. lValtonen et alj200 9) with other aspects receiving less atten- 
tion. One interesting constraint to the models could be obtained 
by monitoring the emission line variability during the outbursts. 
The luminosity of the emission lines could change in response 
to varying continuum radiation or, if the size of the secondary 
black hole orbit and the characteristic broad line region (BLR) 
size are of the same order, the secondary black hole could affect 
the velocity field of the BLR. Furthermore, if the BLR velocity 
dispersion was measured, it could be used to constrain the mass 
of the primary black hole in OJ 287. However, there are no pub- 
lished results of spectroscopic monitoring of OJ 287 in the liter- 
ature. Thus we started a project in 2005 at the VLT to monitor 



* Based on observations collected at the European Organization for 
Astronomical Research in the Southern Hemisphere, Chile, propos- 
als 075.B-0527, 076.B-0028, 077.B-0015, 078.B-0030, 079.B-0014, 
380.B-0665 and 381.B-0050. 



the broad line region of OJ 287 during the expected outbursts in 
2005-08 when the two black holes were expected to be at closest 
approach by various models. The main aim of this project was 
to see if any changes could be observed in the BLR of OJ 287 
during the expected closest approach of the binary. 

This kind of monitoring project is quite challenging be- 
cause the emission lines are very weak in BL Lacs by defini- 
tion. There are only two reported detections of broad lines in 
OJ 287 in the literature. IStickel et all (Il989h detected the H/3 
line with a FWHM of 56 A with an equivalent width (EW) of 
1.1 A. Sitko & JunkkmnenI (Il985h also detected a weak HjS, 
but the most notable feature in their spectrum is the prominent 
broad (FWHM ~ 120 A) Ha line. It is impossible to know if 
these detections represent typical levels or very high states of the 
lines. Given the weakness of emission lines in BL Lacs in gen- 
eral, the latter alternative seems more probable as non-detections 
are not always reported. Based on this past information the Ha 
line was chosen for the monitoring because it can be easily de- 
tected with modern spectrographs at sufficiently large telescopes 
if the luminosity of t he lin e is close to what was observed by 
iSitko & JunkkarinenI (Il985l) . The results of this monitoring are 
reported in this paper 

Throughout this paper we use the cosmology Hq - 70 km 
S-' Mpc-', Qm = 0.3 and Qa = 0.7. 



2. Observations and data reduction 

The observations were made with the FORS2 instrument 
(lAppenzeller et al ] |1998 >) attached to the VLT-UTl (Antu) in 
service mode during seven epochs in 2005-08. Throughout the 
whole period the same instrumental setup was used, employ- 
ing the 1028z grism with the OG590 order separation filter The 
wavelength range covered by the 1028z grism is AA 7730-9480 
A and the resolution with the 2"slit was 8.1 A. The observations 
are summarized in Table [T] The exposure time in Table [T] gives 
the total exposure time of 2-5 exposures of OJ 287 made a each 
epoch. 
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Table 1. Summary of the observations. The last column gives 
the signal-to-noise at 8700 A. 



Epoch 


Date 


texp 


FWHM 


Airmass 


S/N 






(s) 


(arcsec) 






1 


2005-04-05 


1200 


0.9 


1.46 


280 


2 


2005-11-19 


1800 


0.9 


1.61 


330 


3 


2006-04-02 


1500 


0.6 


1.41 


400 


4 


2006-12-06 


1800 


0.9 


1.50 


315 


5 


2007-04-04 


2500 


0.7 


1.43 


440 


6 


2008-01-05 


2500 


0.8 


1.41 


560 


7 


2008-04-08 


2640 


1.1 


1.42 


410 



The reduction of the spectra was made with IRAF. First 
a bias frame was subtracted from the two-dimensional CCD 
spectra after which they were divided by a flat field frame 
derived from continuum lamp exposures. A two-dimensional 
wavelength calibration solution was derived from strips of data 
separated by 5 arcsec in the spatial direction using the same 14 
night-sky emission lines in each strip. This solution was used 
to rectify the image before subtracting the night-sky emission 
lines. After the night-sky subtraction the one-dimensional spec- 
trum of OJ 287 was extracted by tracing the spectrum of OJ 287 
on the CCD and integrating the light within a 7.5 arcsec aperture. 
Optimal extraction was used at this phase to obtain maximum 
signal-to-noise. 

The spectra were calibrated using five spectrophotometric 
calibration stai's (LTT 4816, EG 274, LTT 7379, LTT 6248 and 
GD 108) obtained with the same instrumental setup as here dur- 
ing the observations or from the ESO archive observations made 
in 2005-06. The average sensitivity curve was derived form these 
five stars and applied to the spectra of OJ 287. An extinction cor- 
rection was made to each spectrum using a standard atmospheric 
extinction curve and the airmass at the time of the observation. 
Some of the nights were affected by thin clouds, whose absorp- 
tion cannot be accounted for by the above calibration procedure. 

The calibrated spectra exhibit absorption bands from at- 
mospheric water vapor, which are particularly prominent at 
/1/18 100-8400 A and ^^8800-9300 A and whose strength varied 
from one epoch to another. To correct for the atmospheric ab- 
sorption bands a template was derived from the spectrum of the 
nearby star 10 (iFiorucci & Tostiiri996!) observed simultaneously 
with OJ 287 on 2008-01-05. This template was scaled appropri- 
ately to obtain the best correction for each epoch. During this 
phase it became evident that due to the variable nature of the at- 
mospheric absorption bands they cannot be completely removed 
at their cores. The wavelength regions most affected by this ef- 
fect were excluded from further analysis (see below). 



3. Analysis and results 

Figure [T| shows an example of a calibrated spectrum of OJ 287 
obtained during the monitoring campaign. The spectrum shows 
two noisy areas at /1/18 100-8400 A and /1/18850-9200 A due to 
the imperfect subtraction of the atmospheric absorption bands. 
The spectra also show two areas populated by weak emission 
hues: atA~ 8570 A is the ^^6548,6583 [Nil] + /16563 Ha group 
and at ^ ~ 8780 A the /1/16716,6731 [SII] hues. Horizontal bars 
in Fig.[T|indicate the wavelength areas representing "pure" con- 
tinuum. These areas, 7750-8100 A, 8400-8460 A, 8680-8750 A 
and 8800-8900 A are used in all subsequent work for continuum 
fitting. 



5.5 |-r 




7800 8000 8200 8400 8600 8800 9000 9200 
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Fig. 1. Spectrum of OJ 287 on 2008-01-05. The dashed Hne is a 
power-law fit to pure continuum parts of the spectrum, indicated 
by horizontal bars. The spectral lines identifiable in the spectrum 
have also been marked. 



The overall shape of the continuum follows a power-law 
very well as indicated by the power-law fit in Fig. [1] In all 
seven spectra the maximum deviation from the power-law fit 
is 1%. However, because the main interest of this work is on 
the line variations of OJ 287, a second order polynomial was 
chosen for the continuum fit as it gives about two times smaller 
residuals than the power-law fit and still retains the smooth and 
monotonous shape of the power-law spectrum. Figure |2] shows 
the continuum-subtracted spectra during the seven epochs of the 
monitoring. We note that the continuum level varied by a factor 
of 2.6 during the monitoring. 

During at least two of the epochs (four and seven) a broad 
spectral feature is seen at the location of the Ha-[NII] lines and 
epochs five and six show a hint of a similar feature. However, as 
is evident from Fig.|2l the continuum level is fluctuating around 
the mean level due to calibration errors etc. Because these fluc- 
tuations could be wrongly interpreted as true features, it is im- 
portant to study the statistical properties of the fluctuations more 
closely. 

To study the significance of the spectral features in Fig.|2]we 
employed the method by Sbarufatti et al. (2005). This method 
consists of dividing the spectrum into spectral bins, comput- 
ing the equivalent width (EW) of each bin and deriving the rms 
fluctuations of the EW through the spectrum. Any EW exceed- 
ing a pre-defined rms limit (e.g. 3cr) is taken as a sign of sig- 
nificant emission or a bsorption line warranting further study. 
ISbarufatti et al.l (l2005b computed the continuum level for each 
bin from two adjacent bins, which works fine for narrow lines, 
but very broad feature s may be missed. Thus the method by 
'Sbarufatti et alj ( l2005h was slightly modified by using the fit- 
ted continuum value instead of the adjacent bins as the contin- 
uum estimate. The rms noise in the EW was computed from 
the pure continuum areas marked with horizontal bars in Fig. 
n\ A bin size of 8 A was used in the computations and as in 
ISbarufatti et al. I (l200l any EW value thrice the rms EW was in- 
terpreted as a significant spectral feature. The two areas where 
the atmospheric line correction produces too much noise (Fig.[T]) 
were excluded from the analysis. 

The results of this analysis are shown in Table|2l Column 3 of 
the table gives the derived 3cr EW limit used to detect significant 
features. Columns 3-7 list whether the five narrow emission lines 



K. Nilsson et al: H-alpha monitoring of OJ 287 in 2005-08 



E 
o 

7 

(0 




Table 3. Results of the line fitting. 
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Fig. 2. Continuum-subtracted spectra of OJ 287 during the mon- 
itoring. Individual spectra have been shifted by 10"'^ erg s ' 
cm"^ A"' in vertical direction for clarity. 



seen in the spectrum were formally detected or not. The last two 
columns give wavelength regions where broad spectral features, 
defined as at least three adjacent bins with significant emission, 
were detected and the last column gives any other significant 
spectral features detected. 

Of the five visually identified narrow emission lines only 
/16583 [Nil] is consistently detected throughout the campaign. 
The rest are detected more sporadically, especially the [SII] lines 
are mostly below the 3cr level. There is a clear indication of 
broad emission on top of the [NII]/Hq' lines, which is most natu- 
rally explained as broad Ha emission. This broad emission com- 
ponent is particularly clear at epoch four. In addition to the lines 
in Table |2] the only detected significant features are two weak 
(3cr) absorption features at 8058 and 8738 A. These lines do not 
correspond to any commonly observed spectral lines at the red- 
shift of OJ 287. 

To study the properties and significance of the emission lines 
further, model fits were made to the observed spectra. The nar- 
row lines were modeled by a Lorentzian profile and the broad 
Ha line with a Gaussian profile. Only significant features were 
fitted and the broad Ha was fitted simultaneously with the nar- 
row lines on top of it. Figure|3]shows the result of the line fitting 
on epoch four. Table [3] lists the values of the fitted parameters 
and their associated errors. 

The error bars were estimated by performing model fits to 
500 simulated spectra for each epoch. The simulated data sets 
contained all significant emission features for that particular 
epoch. The noise in the spectra was simulated by performing 
a Fourier transform to the continuum-subtracted pure continuum 
parts of the observed spectra, randomizing the phases and trans- 
forming back to obtain the simulated noise. In this way the statis- 
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FWHM 
(A) 



A, 
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Hey narrow 
6583 [Nil] 
6583 [Nil] 
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6731 [SII] 
6548 [Nil] 
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Fig. 3. Example of the line fitting (epoch four). The data are 
shown as a continuous line and the model fit as a dashed line. 
The lower part shows the residuals of the fit, shifted downwards 
by 0.05 xlQ-^K 



tical properties of the noise were preserved and the fluctuations 
the simulated spectra resembled closely the observed fluctua- 
tions. Examining the distributions of the simulation results it was 
found that for some parameters (e.g. the FWHM of the narrow 
lines) the distribution was not symmetrical, but skewed towards 
higher values. For this reason asymmetric error bars are given in 
Table [3] The quoted error range includes 67% of the simulated 
values, computed separately above and below the median value. 
The broad Ha was considered significant only if the flux of the 
line could be determined to better than 3cr accuracy, where cr is 
the lower error bar. Table [3] lists only significant detections of 
the broad Ha (epochs four and seven). We also computed the 
3cr upper limits for the broad Ha line by creating simulated data 
sets of ~ 300 simulations, each containing increasing contribu- 
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Table 2. Results of the search for significant spectral features. 













Detected? 






Broad 


Other 


Epocli 


JD - 2450000 


EW limit 


[Nil] 


Ha 


[Nil] 


[SII] 


[SII] 


features 


features 






(A) 


6548 


6563 


6583 


6716 


6731 










(A redshifted) 


8550 


8569 


8595 


8769 


8789 






1 


3466.56759 


0.082 


no 


yes 


yes 


no 


no 




absorption at 8738 A 


2 


3693.83049 


0.092 


no 


no 


yes 


no 


no 






3 


3827.53441 


0.072 


no 


no 


yes 


no 


no 






4 


4075.80737 


0.059 


yes 


yes 


yes 


yes 


yes 


emission at 8538-8642 A 


absorption at 8058A 


5 


4194.55954 


0.051 


yes 


yes 


yes 


no 


yes 


emission at 8546-8586 A 


6 


4470.76465 


0.067 


yes 


yes 


yes 


no 


no 


emission at 8546-8570 A 




7 


4564.53963 


0.061 


yes 


yes 


yes 


no 


no 


emission at 8538-8602 A 
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Fig. 4. Line positions (upper panel) and flux density (lower 
panel) of the broad Ho' and /16583 [Nil] lines. In both panels 
the results for the Ha line are shown with closed symbols (upper 
limits with crosses) and the /16583 [Nil] line with open symbols. 
Horizontal lines represent average values. In the upper panel the 
the average value has been subtracted and the plot for the Ha 
moved 30 A up for clarity. 



tion of the broad Ha and determining the line flux which could 
be detected with 3cr significance. 

The rms scatter in the position of the /16583 [Nil] line is 
0.8 A, consistent with the ~ 1 A error bars derived by the error 
simulations. The average redshift derived from the /16583 [Nil] 
line is 0.3056 + 0.0001, consistent with the earlier r edshift z = 
0.306 (.IVliller et al.. 1971 ISifl'CO & Junkkarinenlll985l) . The eiTor 
in the redshift represents the internal precision of the measure- 
ments without systematic elTects, which are expected to be small 
because the wavelength calibration was derived from sky lines 
exposed simultaneously with the OJ 287 spectrum. 

Figure |4] shows the position and flux density of the /i6583 
[Nil] and Ha lines. The /16583 [Nil] fine is not expected to 
show any variations over the timescales considered here, which 
is also confirmed by our results: no variability larger than the 
error bars is observed. However, the error bars of the flux den- 
sities are quite large. The broad //„ line shows no significant 



changes in position or flux between the two epochs it was de- 
tected. As the upper limits show, the data are consistent with no 
change in the broad Ha line flux throughout the whole moni- 
toring campaign. Thus we find no evidence of changes in the 
BLR in our data. However, there is a big diff 'erence to the line 
flux observed by ISitko & Junkkarin"enl (Il985h . who reported a 
Ha line flux of 2.4 x 10" ''^ erg cm"^ s~', which is a factor of 
~ 10 higher than in 2005-08. The difference remains even af- 
ter subtracting the summed contribution of the [Nil] and narrow 
Ha lines, 12 x 10~^^ erg c m~^ s which were unresolved by 
ISifl'CO & Junkkarinenl(ll985h . 

The weighted average of the broad Ha line width is 120 + 
14 A, w hich corresponds to 4200 ± 500 km/s at the redshift of 
OJ 287. ISitko & .Tunkkarinenldl985h do not give the line width, 
but estimating from their Fig. 2 it is ~ 120 A, i.e. the same as 
observed here. 



4. Discussion 

The most notable result of our monitoring is the dramatic de- 
crease in broad Ha luminosity between 1984 and 2005-08. In 
Fig.|5]we show the continuum and broad Ha line flux of OJ 287 
in 1980-2010. In the upper panel three outburst seasons of 
OJ 287 can be seen: 1982-84, 1994-95 and 2005-08. The high 
Ha flux appeared just after the strong outbursts in 1982-83, 290 
days after the second peak. Unfortunately, the sampling is too 
poor to draw any conclusions on the possible connection be- 
tween continuum and BLR luminosity in OJ 287. Such a con- 
nection would not even be necessarily expected because in BL 
Lac objects the vast majority of the continuum arises from the jet 
and is highly beamed by relativistic eff'ects due to the small angle 
between the line of sight and the jet axis. Thus the observed con- 
tinuum variations are not necessarily connected to the changes of 
the continuum source illuminating the Ha emitting clouds, most 
likely the inner parts of the accretion disk. It is nevertheless an 
intriguing observation that the high line luminosity in 1984 was 
observed right after a luminous continuum outburst. Given the 
weakness of the broad Ha line we are unable to study a possible 
connection between BLR velocity field and the suggested per i- 
astron of the secondary black hole in 2005-07 (IValtonenll2007h . 

The average line luminosities in 2005-08 were 
log(L/ergs-') = 41.8, 41.1, 41.2 and 40.8 for broad Ha, 
narrow Ha, ^6583[NII] and ^6548[NII], respectively. The Ha 
luminosities of Seyfert galaxies and quasars are generally much 
higher as shown by he compariso ns in Table IH and Fig. 6] The 
quasar data were obtained from Schneider et alj (12007 ), who 
give data of 77 429 quasars f rom the Fifth Data Releas e of th e 
Sloan Digital Sky Survey dAdelman-McCarthv etsU l2007h . 
We selected quasars at the redshift interval z - 0.30 - 0.31 
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Fig. 5. Continuum (upper panel) and Ha (lower panel) flux of 
OJ 287 in 1980-2010. In the lower panel crosses with arrows 
mark 3cr upper limits and solid circles detections. The error bar 
of the leftmost point in the lower panel has been arbitrarily set 
to 10%. 

Table 4. Comparison of line luminosities of OJ 287 to those of 
quasars and Seyfert galaxies. 







LogL 




Line 


OJ 287 


Quasars 


Seyferts 


Hcf broad 


41.8 


42.9 ± 0.4 


43.1 ±0.3 


/16548[NII] 


40.8 


41.3 + 0.8 




/16583[NII] 


41.2 


41.4 ±0.7 





(147 quasars) and obtained the Ha line fluxes and widths from 
the SpecLine tables in the SDSS archive. Figure |6] shows only 
quasars with FWHM > 1000 km s ' (101 quasars) . In addition 
we plo t data for 17 z = 0.36 Seyfert galaxies in McGill et alJ 
(l2008h . As can be seen from Fig. |6] the Ha luminosity of OJ 287 
was lower than in typical quasars and Seyfert galaxies by a 
factor of ~ 10 in 2005-08. In December 1984, however, the Ha 
luminosity was comparable to that of quasars and Seyferts. We 
finally note that the luminosities of the two narrow lines /16548 
and /16583 [Nil], however, are comparable to the narrow-line 
luminosities of the quasars in I Schneider et alj (l2007h . 



5. Summary 

We have presented high S/N spectra of the BL Lac object OJ 287 
during seven epochs in 2005-08. Our results can be summarized 
as follows: 

1) We were able to detect five narrow emission lines, 
i^6548,6583[NII], ^6563Ha and ^^6716,6731 [SII] during at 
least one of the epochs and a broad Ha feature during two 
epochs. The luminosities of the [Nil] lines are comparable to 
those in quasars at the same redshift, whereas the broad Ha line 
is a factor of ~ 10 less luminous than the Ha line in quasars and 
Seyfert galaxies. 
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Fig. 6. Comparison of broad Ha luminosities of OJ 287 
(open circles connecte d by a line), z = 0.30-0.31 quasars from 
[Schneider et al. ( 20071 filled circles) and z - 0.36 Seyfert galax- 
ies from McGill et alj 12008, open triangles). 



2) The luminosity of the broad Ha line was a factor of ~ 10 
lower in 2005-08 than in 1984 when it was reported to have been 
detected the last time. 

3) We do not see any significant change in luminosity, posi- 
tion or width if the broad Ha line between the two epochs it was 
detected in 2005-08. The average width is 4200 km s"', consis- 
tent with the width observed in 1984. Due to the weakness of 
the broad emission lines we are unable to study the connection 
between continuum and line luminosity or between the BLR ve- 
locity field and the proposed periastron of the secondary black 
hole in 2005-07. 
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